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Bat coronaviruses (Bt-CoVs) are thought to be the 
precursors of severe acute respiratory syndrome corona- 
Virus. We detected BI-CoVs in 2 bat species from Trini 
dad. Phylogenetic analysis of the RNA-dependent RNA 
polymerase gene and helicase confirmed them as group 1 
coronaviruses. 


Postion itt reset fo potest 
ly emergent pathogens. Because of their abundance, 
wide distribution, and mobility, bats confer a greater risk 
for zoonotic transmission than other animals (1). Bats have 
long been known as the natural hosts for rabies virus and 
other lyssaviruses and were more recently identified as the 
reservoirs for emerging viruses such as Ebola, Hendra, and 
\Nipah virus (reviewed in 1). The search for the animal res- 
ervoir of the severe acute respiratory syndrome coronavirus 
(SARS-CoV) led to extensive surveys of coronaviruses in 
‘wild and domestic animal populations in China, resulting in 
the detection of a wide variety of novel bat coronaviruses 
(BELCOVs) (2-5). The data suggest that the progenitor of 
the SARS-CoV, and all other coronaviruses in other animal 
hosts, originated in bats (4). Recent reports by Dominguez 
ct al, (6) and Gloza-Rausch et al, (7) confirmed the exis- 
tence of Bt-CoVs outside China, in the United States and 
Germany, respectively. Additionally antibodies reactive 
with SARS-CoV have been detected in African bat species 
(8). We report the detection and characterization of CoVs 
in bats from Trinidad, the southernmost island of the Carib 
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bean archipelago, located 9 km (5.5 miles) off the north- 
eastern coast of South America, 


The Study 

‘A total of 114 bats collected from their natural habitats, 
from December 2006 through July 2007 (8 species from 10 
locations; Table) were euthanized after being deeply anes- 
thetized with 2% xylazine and 10% ketamine administered 
subcutaneously (in some cases after sedation with CO. 
‘The bats were then taxonomically classified on the basis of 
morphology, and the carcasses were stored at -70°C until 
used, For sampling, the carcasses were thawed at 4°C for 
3-4 hours, then oropharyngeal and anal samples were taken 
with Dacron-tipped swabs that were then placed in RNAI- 
ater (Ambion, Austin, TX, USA) and stored at ~20°C until 
used. 

CoV detection and sequencing were conducted as pre- 
viously described (4). Briefly, viral RNA was extracted 
from swabs by using the QLAamp viral RNA minikit (QUA- 
GEN, Westburg, the Netherlands) and used as the template 
for reverse transcription-PCR (RT-PCR) detection of the 
CoV RNA-dependent RNA polymerase (RdRp) gene (9). 
Primers based on the RdRp gene, conserved for all known 
coronaviruses, were then used for RT-PCR detection. 
‘The RdRp PCR products were gel purified by using the 
QUAquick PCR purification kit (QUAGEN) and sequenced 
to confirm virus species. RNA from samples positive for 
coronavirus was then used for CDNA synthesis by using 
random hexamer, gene-specific, and oligo(dT) primers 
‘The RdRp gene and 1b open reading frames, including the 
helicase (HEL) domain, were then sequenced, also as pre- 
viously described (4). Sequences derived from this study 
were deposited in GenBank (accession nos. EU769557 and 
EUT69558), 

Sequences were aligned with previously published 
CoV sequences from GenBank by using ClustalX (hitp: 
bips.u-strasbg. r/fr/ Documentation/ClustalX) then manual- 
ly aligned by using the Se-Al program (http:/tree.bio.cd.c. 
uk/software/seal). The GenBank accession numbers of all 
sequences used are noted in the taxon names in Figures 1 
and 2. The RdRp sequences were trimmed to equal length, 
which created 2 datasets of 780 bp (n = 40) and 378 bp (n 
= 45). The latter included Bt-CoV sequences from North 
America and Germany that were too short to be included 
in the first dataset. A third dataset comprised an alignment 
of the HEL domain (n = 46) trimmed to 1,797 bp. Maxi- 
mum likelihood (ML) phylogenies were inferred under a 
General Time Reversible (GTR + I+ 1) model, which was 
identified as the best-fit model of nucleotide substitution 
using MODELTEST version 3.7 (10). Bootstrapping was 
performed to assess the robustness of tree topologies by 
using 1,000 replicate neighbor-joining (NJ) trees under the 
ML substitution model. All analyses were performed with 
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‘Table, Number and location of bat species collected and tested from December 2006 thraugh July 2007, with bat coronavirus species 


status 


Species (amy) 


Tocaiion (ay 


‘Total no, Tested (90- POSTIVE) 


‘Carolla perspiailata (Phylostonvdae) “Arima @), Fyzabad (2), Tabagute Gl) 5 (0) 
Glossophaga soricina (Phylostomidae) Couva (12)*, Tabaquite (7) 21) 
‘Noctlio leporinus (Noctlionidae) ‘Cawa (6) 500) 
Desmodus rotundus (Phyllostomidae) __—Fyeabad (3), Mocne Diablo (3), Rousillac (2), La Brae (7) 140) 
Preronones parneli (Mormoopidae) ‘Tabaquite (29), Wallerfield (2) 31 (0) 
‘Molossus major (Molossidae) Talparo (25) 25 (0) 
‘Mormoops sp. (Mormaopidae) ‘Tamana (1) 1(0) 
Phyllostomus hastatus (Phylastomidae) Waller (11) 110) 
Total ‘Wlocations: 8 species 14 


Incite Grup Wa WHET Co PSS Bas STE 


PAUP® version 4.0b (Sinauer Associates, Inc., Sunderland, 
MA, USA). 

CoV RNA was detected in | of 21 Glossophaga soric- 
ina and 1 of 5 Carollia perspicillata bats tested. The latter, 
designated Bt-CoV/Trinidad/IFY2B, was an adult male 
from the area of Fyzabad, and virus was detected in anal and 
oropharyngeal swabs. The infected G. soricina specimen 
(BL-CoV/Trinidad/1CO7B) was an adult female from the 
area of Couva, and the virus was detected in the anal swab 
only. Subsequent sequencing was performed on virus from 
the anal swabs. A total of 3,905 bp of RdRp from Bt-CoVi 
Tri/1CO7B and 5,160 bp from Bt-CoV/Tri/IFY2B were 
sequenced, in addition to 1,782 bp from the HEL domain 
of both samples (full sequence data not shown). Across the 
shared 3,798-bp region of RdRp, divergence was 28.7% 
at the nucleotide level and 15.5% at the amino acid level, 
Similarly, the HEL domains of the 2 viruses showed 27.3% 


a2) 


24) 


26) 


194) 


divergence at the nucleotide level and 13.3% at the amino 
acid level, 

ML phylogenies inferred from the RdRp genes of 
40 viruses (780 bp) and the HEL domains of 46 virus- 
es (1,797 bp) are shown in Figure 1. In both eases, the 
5 groups proposed by Tang et al. (4) on the basis of NJ 
trees, and the lineage containing the recently reported 
novel CoV sequences (11), were strongly supported with 
bootstrap values >95% in all cases. In each case the Trini- 
dadian sequences clustered with group | CoVs within a 
clude containing all other group 1 bat and human CoVs 
as well as porcine CoV. To determine the phylogenetic 
relationship between the Trinidadian Bt-CoVs and North 
American Bt-CoVs, for which only relatively short se~ 
quences for RdRp were available, we inferred a second 
RdRp ML phylogeny based on a 378-bp fragment (Figure 
2), When this shorter fragment was used, 5 groups were 


204) 


Figure 1. Maximum lkeliiood tees of coronaviruses based on A) 780-bp fragment of the RNA-dependent RNA polymerase gene and 8) 
1,797 bp of the helicase (HEL) domain of open reading frame 1b. Trees were inferred under the General Time Reversible (GTR +, + 
1) model by using PAUP* version 4.0b (Sinauer Associates, Inc., Sunderland, MA, USA). Bootstrap support values >80% are Indicated 
Previously defined phylogenetic groups and a putative navel group (20) are delineated by the bars on the right. The numbering of these 
‘groups is as described in the eighth report of the Intemational Committee on Taxonomy of Viruses with the alternative grouping propased 
by Tang etal. (4) in brackets. Trinidadian bat coronavirus sequences are highighted in red. GenBank accession numbers are noted in 
parentheses. Scale bars indicate number of nucleotide substitutions per ste. 


Emerging Infectious Diseases + www cde.govled + Vol. 14, No. 12, December 2008 p91 


DISPATCHES 


wo 


Figure 2. Maximum lkelhood tree of coronaviruses based on 378- 
bp fragment of the RNA-dependent RNA polymerase gene. The 
‘wee was inferred under the General Time Reversible (GTR ++!) 
bby using PAUP* version 4.00 (Sinauer Associates, Inc, Sunderland, 
MA, USA). Trinidadian bat coronavirus (Bt-CoV) sequences are 
hightghted in red and North American BI-CaV in blue. Previously 
defined phylogenetic groups and a putative navel group (20) are 
dolineated by the bars onthe right. The numbering of these groups 
|g as described in the eighth report of the international Commitee 
‘on Taxonomy of Viruses with the alternative grouping proposed by 
‘Tang et al. (4 in brackets. Bootstrap suppor values for groups 1a, 
4b, 2a-c, 3, and the lineage containing Trinidadian 8t-CoVs are 
‘shown. GenBank accession numbers are noted in parentheses. 
‘Scale bars indicate number of nucleotide substiutions per sit 


again defined, but bootstrap support for the putative group 
5 (4) was lower. The level of divergence between Trinida- 
dian sequences was notably higher than among the North 
American (6) and German sequences (7). 


Conclusions 

Our detection of RNA from group 1 CoVs in Trini- 
dadian bats shows that Bt-CoVs have a wider distribution 
than previously suspected and is added support for bats as 
the original host species for these viruses. Group 1 CoV 
form 2 well-supported clades designated 1a and 1b (12). 
The Trinidadian bt-CoV clustered within the latter clade, 
\hich contains all other group 1 BL-CoVs, including those 
from Germany and North America, and the 3 known group 
1 human CoV associated with respiratory illness (13-15). 
Despite the geographic proximity of the bats from which 
the Trinidadian Bt-CoV sequences were derived—Couva 
and Fyzabad are 28 km (17 miles) apart, and Trinidad has 
‘an area of only 4,769 km? (1,864 square miles)—they are 
relatively highly divergent. This divergence might reflect 
virus adaptations to different host species; however, more 
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data would be needed to confirm this. Given the mobility 
of bats, the possibilty of the viruses having different geo- 
graphic origins (perhaps even from outside Trinidad) can- 
not be ruled out. Further work on CoV diversity in Trinidad 
and the rest of the Americas, as well as on the ecology and 
bchavior of susceptible bat species, is needed to understand 
the origins, evolution, and dispersal of these viruses. 
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